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This  work  is  a  study  of  the  Stark  broadening  of  spectral  lines 
emitted  by  highly  charged  ions  in  a  hot,  dense  plasma.  The  line 
broadening  calculations  of  Tighe  and  Hooper,  in  which  the  dipole 
approximation  was  used  for  the  interaction  between  the  radiating  ions 
and  perturbing  electrons,  are  extended  by  retaining  the  full  Coulomb 
radiator-perturbing  electron  interaction  .  Electron  broadening  is 
treated  to  second  order  in  the  radiator-perturbing  electron  interaction; 
ion  broadening  is  treated  by  a  static  ion  microfield  probability 
distribution  .  Perturbing  electrons  are  treated  quantum  mechanically 
through  Coulomb  wavefunctions  to  account  for  the  charged  radiator . 

Line  profiles  from  this  calculation  are  slightly  broader  than  those 
calculated  using  the  dipole  approximation  .  The  near  agreement  of  these 
results  is  partly  fortuitous  in  that  the  dipole  approximation 
overestimates  the  dipole  part  of  the  Coulomb  interaction  for  perturbers 
close  to  the  radiator,  partially  compensating  for  the  neglected 


multipoles  of  the  interaction .  This  full  Coulomb  calculation  of  the 
dynamic  line  shift  due  to  perturbing  electrons  results  in  a 
significantly  smaller  shift  than  that  obtained  when  the  dipole 
approximation  for  the  radiator-perturbing  electron  interaction  is 
used .  For  an  Argon  Lyman  fj  line  calculated  for  the  plasma  conditions  of 
10   electrons/cm  and  800  eV,  this  shift  is  about  0.03  Ryd  toward  lower 
energies,  and  it  causes  no  noticeable  asymmetry  in  the  line.  The  effect 
of  the  dynamic  shift  on  the  line  is  considerably  smaller  than  that  of 
the  static  (plasma  polarization)  shift  calculated  by  Skupsky,  which 
causes  an  asymmetry  in  the  p  line,  with  the  red  wing  having  the  greater 
intensity.  For  an  Argon  Lyman  |3  line,  the  significance  of  ion-radiator 
and  electron-radiator  perturbation  matrix  elements  between  the  initial 
radiator  states,  principal  quantum  number  =  3,  and  the  states  of  the 
nearest  adjacent  energy  level,  principal  quantum  number  =  4,  is 
examined  .  The  electron-radiator  matrix  elements  are  found  to  give  a 
small  additional  broadening  of  the  line;  the  ion-radiator  matrix 
elements  are  found  to  cause  a  significant  nonlinear  Stark  shift  of  the 
radiator  energy  levels.  This  shift  causes  the  line  to  be  asymmetric, 
with  the  blue  peak  having  greater  intensity  than  the  red  peak  and  the 
red  wing  having  greater  intensity  than  the  blue  wing  . 


CHAPTER  I 
PLASMA  SPECTRAL  LINE  BROADENING 


1.1  Introduction 


Plasma  broadened  spectral  lines  have  been  used  for  many  years  in 
determining  the  densities  and  temperatures  of  laboratory  and  astro- 
physical  plasmas  (Griem  1964,  1974;  Baranger  1962;  Cooper  1966;  Smith, 
Cooper,  and  Vidal  1969,  1970).   Recently,  plasma  broadened  X-ray  spectra 
from  highly  ionized  high-Z  elements  (e.g.  neon  or  argon)  have  been  used 
to  determine  the  densities  of  laser  inertial  confinement  plasmas 
(Yaakobi  et  al .  1977,  1979,  1980;  Kilkenny  et  al .  1980;  Apruzese  et  al . 
1981). 

The  purpose  of  this  work  is  to  extend  theoretical  line  shape  calcu- 
lations by  Tighe  (1977)  and  Tighe  and  Hooper  (1976,  1978),  which  were 
used  in  the  studies  of  Yaakobi  et  al .   The  work  of  Tighe  and  Hooper  is 
an  adaptation  of  the  relaxation  theory  of  line  broadening  (Smith  1966; 
Smith  and  Hooper  1967;  O'Brien  1970;  O'Brien  and  Hooper  1974)  to  hydro- 
genic  radiators  of  arbitrary  charge,  appropriate  for  the  high-Z 
radiators  of  laser  implosion  experiments.   In  the  Tighe-Hooper  calcula- 
tions the  dipole  approximation  is  made  for  the  radiator-perturbing 
electron  interaction,  an  approximation  which  is  of  questionable  validity 
at  the  high  densities  obtained  in  some  recent  experiments.   In  this  line 
shape  calculation  we  will  not  make  the  dipole  approximation  but  will 
retain  the  full  Coulomb  radiator-perturbing  electron  interaction. 
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In  this  chapter  we  outline  a  line  broadening  theory  which  is 
obtained  by  use  of  the  kinetic  theory  of  time  correlation  functions 
(Hussey  1974;  Hussey  et  al .  1975),  then  apply  several  approximations  to 
obtain  a  theory  in  which  the  line  width  and  shift  operator  is  expanded 
to  second  order  in  the  radiator-perturbing  electron  potential .  In 
Chapter  II  we  develop  a  computational  form  for  the  full  Coulomb  second 
order  theory .  In  Chapter  III  we  present  Lyman  series  line  shapes 
calculated  numerically  from  the  formalism  of  Chapter  II  .  We  compare 
profiles  calculated  from  this  full  Coulomb  theory  with 
the  results  of  Tighe  and  Hooper,  and  also  compare  this  calculation 
to  a  recent  impact  theory  of  line  broadening  (Griem  et  al  .  1979). 
Calculated  line  shapes  which  include  dynamic  shift  and  inelastic  effects 
are  also  presented  . 

I  .2  The  Line  Shape 
The  power  spectrum  emitted  by  one  type  of  ion  in  a  plasma  can  be 
written  in  terms  of  the  ensemble-averaged  radiation  emitted  by  one  ion 
of  that  type  (the  radiator)  as  (Griem  1974) 
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P(<A>=^3  I      |<a|J|b>|2p  6(u>-w )  (1.2.1) 

3c   ab 


where  p  is  the  probability  that  the  radiator-plasma  system  is  in  the 

.    ■+ 
state  |a>,  d  is  the  dipole  moment  of  the  radiator,  and  u  .  =  (E  -E,  )/fi 

ab     a  b 

The  line  shape  function  I(co)  is  defined  in  terms  of  the  power  spectrum 
by  the  equation 


,  4 
P(u)  =  -^\  I(to), 
3c 


where 


I(o>)  =  I    |<a|d|b>|2  p  5(w-u>  ).  (1.2.2) 

ab 


We  can  express  the  Dirac  delta  function  in  I(co)  as  an  integral  to  obtain 


I(u,)=-k'-~  I      Md|b>|2pa  •1(w-w«b)tdt.         (1.2.3) 
ab 


Here,  the  integrand  for  negative  values  of  t  is  equal  to  the  complex 
conjugate  of  the  integrand  for  positive  values  of  t,  so  we  can  rewrite 
this  equation  in  terms  of  an  integral  from  zero  to  infinity: 

I(o))  -  -±  Re  J"     I      |<a|d|b>|2  p  ei(u~uab),:dt        (1.2.4) 


if    \        l   v,      f°°  V   iwt/uiti  x  y    i   -iHt/ti  -t   iHt/-fi,,v,^ 
I((j)  =■  —  Re  J    I     e   <b|d|a>«<a|  pe       de     |b>dt  , 

ab 

(1.2.5) 


where  H  is  the  Hamiltonian  for  the  radiator-plasma  system  and  p  is  the 
canonical  density  operator,  p  =  e   /Tr  e    .   This  can  be  written  as 


I(w)  =  —  Re  Tr  J   dt  e    d«p  e       d  e   '         (1.2.6) 

ii       •'o 

where  Tr  represents  a  trace  over  states  of  the  radiator-plasma  system. 

In  Eq.  (1.2.1)  we  have  neglected  Doppler  broadening  of  the  line 
shape  due  to  the  motion  of  the  radiator.   This  will  be  included 
approximately  at  the  end  of  the  calculation  by  convoluting  the  line 
shape  with  a  Doppler  profile  based  on  a  Maxwell  velocity  distribution. 
This  is  an  approximate  treatment  of  Doppler  broadening  based  on  the 
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assumption  that  the  change  of  momentum  of  the  radiator  is  negligible 
during  the  time  of  radiation  (Smith  et  al .  1971a,  1971b).   The  Doppler 
profile  is  given  by 

-(to-co  )  /y 
ID(w)  =— J— e     '      ,  (1.2.7) 


where 


2   2kT   2 


Y   =  — r  u. 


Mc 


2  o  ' 


M  is  the  mass  of  the  radiator,  and  u)  is  the  frequency  of  the 
unperturbed  transition. 

1.3  Time  Scales 


The  line  shape  is  given  by  the  Laplace  transform  of  the  radiator 

dipole  autocorrelation  function  (Eq.  1.2.6).   Since  the  transform  has 

the  property  AcoAt  <  1 ,  where  Ato  is  the  frequency  separation  from  the 

line  center,  the  time  dependence  of  a  perturbation  of  the  radiator  which 

changes  significantly  in  time  t  will  be  evidenced  by  the  the  part  of  the 

line  for  which  Aw  <  1/t.   For  Aw  »  1/t,  the  perturber  can  be  considered 

as  static.   This  part  of  the  line  corresponds  to  radiation  from  initial 

radiator  states  having  average  lifetimes  much  shorter  than  t.   In  the 

plasma  line  broadening  problem  there  are  two  characteristic  perturbation 

times,  x.  ~  1/to  .  for  perturbations  due  to  ions  and  t  ~  l/co   for 
l      pi  e      pe 

perturbations  due  to  electrons,  where  to   and  to  .  are  the  electron  and 

pe      pi 

ion  plasma  frequencies. 

For  the  plasma  conditions  which  we  will  examine,  the  half-widths  of 

the  first  few  Lyman  lines  fall  within  the  region  Ato  <  to   ,  so  we  must 

pe 

treat  the  perturbation  due  to  the  electrons  as  a  dynamic  process.   In 
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the  case  of  the  ions,  the  region  Aco  <  u  .  corresponds  to  only  a  small 
part  of  the  line  center.   We  will  therefore  ignore  ion  dynamic  effects 
and  make  the  approximation  that  the  perturbing  ions  and  radiator  are 
static  and  that  the  ions  perturb  the  radiator  through  a  static  electric 
field.   This  approximation,  which  is  used  in  many  line  broadening 
calculations,  is  known  as  the  Static  Ion  Approximation.   Lyman  line 
shapes  which  include  ion  motion  effects  have  been  calculated  for 
hydrogen  plasmas  (Greene  1979;  Seidel  1980).   There,  the  inclusion  of 
ion  motion  was  found  to  cause  additional  broadening  near  the  line 
center.   The  peak  of  the  Lyman  alpha  line  was  significantly  lowered,  the 
peaks  of  the  Lyman  beta  line  were  lowered  slightly  and  the  central  dip 
was  raised. 

1.4  The  Ion  Microfield 
In  this  section  we  apply  two  approximations  which  will  permit  us  to 
write  the  line  shape  expression  in  the  form 

I(u))  =  /"  P(e)J(u,e)  (1.4.1) 

where  P(e)  is  the  low  frequency  ion  microfield  probability  (Baranger  and 
Mozer  1959;  Mozer  and  Baranger  1960;  Hooper  1968;  O'Brien  and  Hooper 
1974;  Tighe  and  Hooper  1977)  and  J(u),e)  is  the  line  shape  due  to  the 
electron  broadening  of  radiators  in  the  presence  of  the  ion 
microfield  e.   The  Hamiltonian  for  the  radiator  and  plasma  is 


H=H+H.  +H+V.+V   +  V.  (1.4.2) 

r    i    e    ri    re    xe 
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where  1^,  E±,   and  Hg  are  the  Hamiltonians  for  the  radiator,  ion,  and 
electron  systems,  Vri ,  Vre ,  and  Vi   are  the  potential  energies  between 
the  radiator  and  ions ,  radiator  and  electrons ,  and  ions  and  electrons  . 
The  effect  of  the  interaction  Vie  is  to  correlate  the  ions  and  electrons 
such  that  the  average  electron  density  is  greatest  near  the  ions  . 
Therefore ,  the  field  at  the  radiator  due  to  the  electrons  has  a  low 
frequency  component  which  tends  to  partially  cancel  the  ion  field  there  . 
If  the  average  electron  spacing  is  much  less  than  the  Debye  length,  i.e. 


-1/3       _  ,        kT   a/2  ..        ,. 

ne     «   \  ~   I  ~ 2  j    »  (1.4.3) 

4im  e 


where  n  is  the  electron  density,  and  quantum  mechanical  effects  are  not 
significant  (Iglesias  1982)  ,  we  can  treat  this  correlation  effect  to  a 
good  approximation  by  considering  the  low  frequency  microfield  at  the 
radiator  to  be  due  to  Debye  screened  ion  fields  (Baranger  1962)  . 
Similarly,  we  can  consider  the  ion-ion  interaction  to  be  screened  by  the 
electrons.  With  the  effects  of  V^e  approximated  by  screening,  we 
eliminate  it  from  the  Hamiltonian  to  obtain 


H  =  (K.  +  VS  +(().)  +  (H  +  H  +  V   +  e  e  «x  ) 
l    x    ri      r    e    re      i   r 


=  H!  +  H'  ,       .  (1.4.4) 


where  H'  is  the  Hamiltonian  for  the  ions  and  H'  is  the  Hamiltonian  for 
the  radiator  and  electrons.  Here,  K.  is  the  ion  kinetic  energy,  V.   is 
the  shielded  ion-ion  potential  energy,  H  is  the  Hamiltonian  for  an 
electron  gas ,  <J>  ,  is  the  monopole  part  of  the  shielded  ion-radiator 
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potential,  e.  is  the  shielded  ion  field, 
1 

Z.ex.       x.    -x./X._ 
J    x.  D 

and  x  and  x.  are  the  positions  of  the  radiator  electron  and  the  i'th 
ion  with  respect  to  the  radiator  nucleus,  which  we  choose  as  the  origin 
of  our  coordinate  system.   We  will  neglect  terms  of  higher  order  than 
dipole  in  the  radiator-ion  interaction.  This  approximation  is  much 
better  for  the  ions  than  for  the  electrons  in  the  case  of  high  Z  per- 
turbers  since  the  ion  density  is  1/Zj_   times  the  electron  density  and 
since  the  large  ion-radiator  repulsion  tends  to  keep  the  ions  and 
radiator  apart.   However,  for  deuterium-tritium  plasmas  with  a  small 
amount  of  high-Z  impurity,  terms  of  higher  order  than  dipole  in  V  ^  may 
have  a  significant  effect  on  the  line  shape. 

We  also  use  the  static  ion  approximation.   We  assume  that  the  ions 
do  not  move  significantly  during  the  lifetime  of  the  initial  radiator 
state,  so  we  make  the  approximation 


am 

Hi  ■?-=■  =  [H!,h1  =  0.  (1.4.6) 

ot    L  1   J 


Then  the   time   evolution   and   density  operators   can  be   factored, 
iHt/fi  iH't/fi     iH't/fi 


(1.4.7) 


and 

P   =  P±Pre    ,  (1.4.8) 

where 

-pH'  -pH'. 

p±  =  e        1/Tr±   e        X  (1.4.9) 


and 


pre   =  e-PH'/Trre   e"PH'     .  (1.4.10) 


The   line   shape   then   is   given  by 

t(    \        1   d     t  C°  j-      i(jJt   1  t.  -iH't/fi  t     IH't/tt 

I(cj)    =  —  Re  Tr .  p .    I      dt   e  d»Tr        p       e  d  e  . 

te  i  i  '  o  re     re 

(1.4.11) 

Next  we  multiply  this  equation  by  6(e-e  )  and  integrate  over  e  to  obtain 

I(a>)  =  /  Q(e)J(w,£)de  , 
where 

Q(e)  =  Trip16(e-ei)  (1.4.12) 

is  the  probability  of  finding  the  ion  microfield  e  at  the  radiator, 

r,        N   1  D  _    f»  ..   iut  *      -iH'(e)t/tt  t  iH'(£)t/h 
J(w,e)  =  —  Re  Tr   J   dt  e    d«p   e  d  e   v  '     , 

n      re  Jo  rre  ' 

and 

H'(e)  =  H  +  H  +  V   +  e£«x   . 
r    e    er       r 

Since  we  consider  an  isotropic  plasma,  we  can  define  the  probability  of 
finding  the  ion  microfield  of  -magnitude  e  at  the  radiator  by 
P(e)  =  4n;£2Q(e)  .  Then 

I(oj)   =  J"  P(e)J(o),e)de  (1.4.13) 

and 


1  o_  ^_        r"  ^   ~iwt  t..       _-iH'(e)t/fi  t  _iH'(e)t/-n 

(1  .4.14) 


t/        \        1   n     ^  f°°    ,        iut  }  -iH'(e)t/'h  -> 

J(co,e)    =  —  Re   Tr        J      dt  e  d«p       e  d 

ti  re   J  o  're 


where 


H»(e)  =  H  +  H  +  V   +  eez 
r    e    er      r 


=  H(r)  +  H  +  V   . 
e    er 


For  convenience,  we  have  chosen  the  z-axis  to  be  in  the  direction  of  the 
field  e. 

I  .5  Kinetic  Theory 

Here  we  follow  the  work  of  Hussey  (1974)  and  Hussey  et  al  .  (1975) 
to  obtain  an  expression  for  J(oj,e)  in  terras  of  a  width  and  shift 
operator,  M(w)  .  The  advantage  of  this  form  is  that  results  obtained  by 
approximate  treatments  of  M(u)  are  much  better  than  those  obtained  when 
similar  approximations  are  applied  directly  to  J(w,s)  . 

We  define  the  Liouville  operator  L  in  terms  of  a  commutator  of  the 
Hamiltonian  H'(e), 

Lf  =  [H\f]/fi  (1.5.1) 

and  write  J(co,e)  as 


J(u,e)  =-Re  Tr-   feiut  o*.p   e"iLt  d~  .         (1.5.2) 
%  re  ;o        Kre  v     ' 


We  then  define  the  radiator  and  s-perturber  functions 


nS  F(r,l,...,s;t)  .  -^  Tr^  #  #  ^   p^  e"iLt  I  ,   (1.5.3) 
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where  Trs+j_..#N  is  a  trace  over  perturbing  electrons  s+1  through  N  and  n 
is  the  electron  density.   We  can  then  write  J(w,£)  in  terms  of  the  s=0 
function  F(r;t)  or  its  Laplace  transform  F(r;(jj)  as 


J(u,e)  =  -  Re  Tr  /"  dt  eiut  d.f(r;t) 
%  r  J  o  '  ' 


=  -  Re  Tr   3'F(r;w)  .  (1.5.4) 


If  we  take  a  partial  time  derivative  of  Eq.  (1.5.3)  we  obtain 


[|j-+  iL(r,l,...,s;t)]f(r,l»...,s;t)  = 


-  in  Tr   [LjCr.s+l)  +  £  L  (i,s+l )]f(r , 1 , . . . ,s+l ;t)  , 
i=l 

(1.5.5) 

where 


L(r,l,...»s;t)f  =  [H(r)  +  He(l,...,s)  +  Ver( 1 , . . . ,s) ,f ]  , 
LjCr.Df  =  [VjCr.D.f]/*  , 

L2(l,2)f  =  [V2(l,2),f]/fi  , 

V,(r,l)  is  the  interaction  between  the  radiator  and  perturber  1,  and 
V"2(l,2)  is  the  interaction  between  perturbers  1  and  2.   Equation  (1.5.5) 
relates  the  s-perturber  function  F(r , 1 , . . . ,s ; t)  to  the  s+1-perturber 
function  F(r , 1 , . . . ,s+l ; t) ,  so  we  have  an  infinite  hierarchy  of  equa- 
tions, the  BBGKY  hierarchy  (Bogoliubov  1962).   We  will  formally  close 
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the  hierarchy  to  obtain  an  equation  for  F(r;co),  and  therefore  J(to,e),  in 

terms  of  a  width  and  shift  operator,  M(w),  to  be  defined  later. 

We  define  the  radiator  and  s-perturber  reduced  distribution 

functions , 

nSf(r,l,...,s)  =  ,..N'    Tr  ..    „p    ,  (1.5.6) 

(N-s)!    s+1  . .  .N  're 

and  the  time  evolution  operator 


U(r'1'""S;t)  =^r-T7Tr^"-NPree"iLtf(r'1"-"S)"1 
n  (N-s)! 


(1.5.7) 


such  that 

F(r,l,...,s;t)  =  U(r,l,...,s;t)F(r,l,  ...,s;t=0) 

=  U(r,l,...,s;t)f(r,l,...,s)d\         (1.5.8) 

The  closure  relation  is  obtained  by  eliminating  d  from  the  s=0  and 
arbitrary-s  members  of  the  Laplace  transform  of  Eq .  (1.5.8), 

F(r,l,...,s;a>)  =  U(r  ,1 ,  ..  .,s  ;w)f  (r  ,1 ,  . .  .,s)f  (r)_1U(r  ;w)_1F(r  ;w)  , 

(1.5.9) 
where 


U(r,l,...s;w)  --^ Trg+1   N  prei(w-L)  *f  (r  ,  1  ,  . .  .,s)  l. 

n  (N-s)! 


(1.5.10) 
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This  equation  for  s=l  is  substituted  for  F(r,l;u)  in  the  Laplace 
transformed  first  equation  of  the  hierarchy  to  obtain 


[-iu  +  iL(r)]F(r;cj)  -  f(r)<T  = 

-inTr1L1(r,l)[u(r,l;co)f(r,l)f(r)"1U(r;(i))~1]F(r;(i)).      (1.5.11) 

Then  we  separate  the  term  in  brackets  on  the  right  side  of  Eq.  (1.5.11) 
into  frequency  dependent  and  frequency  independent  parts  by  defining 


U(r,l,...,s;u)f(r,l,...,s)f(r)  Vrju)  l 


K(r,l s;w)  +  f (r ,  1 , . . .  ,s)f (r)  l  (1.5.12) 


and  solve  for  F(r;u>)  to  obtain 


f(r;co)  =  i[cj  -  L(r')  -  B  -  M(o))]"1  f(r)l  ,  (1.5.13) 


where 


B  =  nTr1L1(r,l)f(r,l)f(r)  l    ,  (1.5.14) 


M(u)  =  nTr1L1(r,l)K(r,l;w)  .  (1.5.15) 


Similarly,  we  can  use  the  second  equation  of  the  hierarchy  (Eq.  1.5.5) 
to  obtain  for  K(r,l;w)  (Hussey  1974;  Hussey  et  al.  1975): 
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l-lr/_  -S  *,         ,**,  s-1 


K(r,l;u)  =  [u-L(r.l)  -  V(r  ,1  ;co)  ]_if  (r  ,1)  £  (r  ,l)f  (r)"1    (1.5.16) 


where 


-1 


and 


V(r,l;u)  =  -nf(r,l)f(r)  1  Tv^^r  t2)?n 


+  nTr^L^r^)  +  L2(l,2)]G(r,l,2;u)G(r,l;to) 


3C(r,l)  =  LjU.l)  +  nf(r,l)-1  Tr2[f(r,l,2) 


-  f(r,l)f(r)~1f(r,2)]L1(r,2)  .  (1.5.17) 


The  operator  P12  acts  to  the  right  to  change  the  argument  1  to  2.   We 
use  Eq.  (1.5.16)  in  Eq.  (1.5.15)  to  obtain 

M(u)  =  nTr1L1(r1l)[u-L(r,l)  -  V(r ,1 ;w)]f (r ,1)3? (r ,l)f (r)"1  . 

(1.5.18) 

With  Eqs.  (1.5.4)  and  (1.5.13)  we  can  now  write  J(w,e)  as 

J(w,e)  =  -  I  Im  Tr  ■  d*  [w-L(r)-B-M(u>)  ]_1  f(r)d  .       (1.5.19) 

Here,  the  effects  of  the  plasma  electrons  on  the  line  shape  are 
contained  in  B  and  M(w).   The  operator  3  is  a  mean  field  or  Hartree- 
Fock  type  term  which  represents  the  shift  of  radiator  energy  levels  due 
to  polarization  of  the  plasma  by  the  radiator.   The  real  and  imaginary 
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parts  of  the  operator  M(oj)  represent  the  dynamic  shift  and  width, 
respectively,  due  to  finite  time  collisional  effects.   The  operator 
j£(r,l)  is  the  interaction  of  the  radiator  and  perturber  1  statically 
screened  by  the  other  perturbing  electrons. 

1.6  Second  Order  Theory 

In  this  section  we  restrict  the  width  and  shift  operator,  M(w)  ,  to 
second  order  in  the  radiator-perturbing  electron  interaction,  V,(r,l). 
We  also  make  use  of  other  approximations  which  are  common  to  many  line 
shape  calculations:  the  neglect  of  interactions  between  perturbing 
electrons  (The  correlation  effects  thus  neglected  will  be  reintroduced 
in  an  approximate  manner  later  in  the  calculation),  the  No-Quenching 
Approximation  (NQA) ,  and  the  No  Lower  State  Broadening  Approximation 
(NLBA). 

The  neglect  of  perturbing  electron  interactions  permits  us  to  write 
the  electron  Hamiltonian  as  a  sum  of  single-particle  Hamiltonians, 


He  =  I      H(i)  .  (1.6.1) 

i 


Since  the  radiator  is  at  the  origin  of  our  coordinate  system,  the  long 

2 
range  monopole  part  of  the  radiator-perturber  interaction,  -(Z-l)e  /x  , 

is  a  function  of  only  the  perturber  coordinate  x  .   Therefore  we  include 

it  in  the  one-electron  Hamiltonian,  H( 1 ) ,  by  redefining  H(l)  and  V,(r,l) 

as 

2 

(1.6.2) 
zm        x 

and 

2      2 

(1.6.3) 


Pi 

H(1)    -2=-     " 

(Z-l)e2 
Xl 

2 

\r     (  r      l\      -               6 

2 
e 

1                                1+        *      1 
X     "X, 

r      1 

Xl 

15 
We  can  now  expand  M(a>)  in  this  short  ranged  V  (r,l)  while  retaining  to 
all  orders  the  effect  of  the  radiator  raonopole  on  the  perturbing 
electrons . 

For  the  plasma  conditions  and  Lyman  lines  which  we  will  consider,  a 
large  part  of  each  line  satisfies  Aw  <  w   ,  where  Aw  is  the  frequency 
separation  from  line  center  and  co   is  the  electron  plasma  frequency. 
The  electron  broadening  of  this  part  of  the  line  is  primarily  due  to 
weak  electron  collisions,  so  we  can  well  approximate  M(oo)  by  retaining 
only  second  order  terras  in  an  expansion  in  V,(r,l)  (Smith,  Cooper,  and 
Vidal  1969).   With  these  approximations  we  have 

M(w)  =  nTTjLjCr.Dtu  -  L(r)  -  L(  l)]"1L1(r  ,l)f  (  1)  ,      (1.6.4) 


where 


and 


f(r)  -  .-P^r),    -pH(r)  Q^5) 


f(l)  -f.-^/Tr^*1'  .  d.6.6) 

In  Appendix  A  we  discuss  the  NQA  and  NLBA,  and  show  that  as  a  conse- 
quence of  these  approximations,  M(o>)  has  the  matrix  form  (Tighe  1977) 


M(w).it  = 


-i^Tr,  fat  I   eiAwV...e-iH(1)t/V...,eiH(1)t/*f(l)  , 
,2   1  'o   h„  n  ii*  K    J    ' 

h  i 

(1.6.7) 
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where  Aco  =  w-(E  ,-E  )/Ti   and  the  subscript  i  (f)  represents  initial 
(final)  states  for  the  particular  Lynan  line  to  be  calculated.   We  will 
partially  include  the  effects  of  quenching  later  in  this  work. 

This  form  of  M(u))  with  the  dipole  approximation  for  V,  (r,l)  was 
used  by  O'Brien  (1970)  and  Tighe  (1977)  in  calculating  Lyman  line  shapes 
for  charged  radiators.   Here,  we  will  retain  the  full  Coulomb 
interaction,  V,(r,l),  in  calculating  Lyman  series  lines. 


CHAPTER  II 
THE  FULL  COULOMB  RADIATOR-PERTURBER  INTERACTION 

II. 1  The  Line  Width  and  Shift  Operator 
The  purpose  of  this  chapter  is  to  develop  a  computational  form  for 
M(co)  in  which  we  make  no  approximations  for  the  radiator-perturbing 
electron  interaction,  Eq.  (1.6.3).   We  will  calculate  matrix  elements 
of  M(o))  in  the  spherical  representation,  |nlm>,  where  n  is  the  principal 
quantum  number  of  the  initial  level  of  a  Lyman  transition.   The  perturb- 
ing electron  Hamiltonian,  H(l),  is  that  of  a  particle  in  an  attractive 
Coulomb  potential,  so  we  will  use  the  Coulomb  wavef unctions  to  evaluate 
the  trace  in  M(<jj)  (O'Brien  1970).   The  eigenfunctions  of  H(l)  are  (Alder 
et  al.  1956) 


where 


<tlh  '   L  77T372  S   e10«'k)Yi>  V^^-'Yl.^ 

xm  ( An) 


aU,k)  =  argrU  +  1  +  in), 


_  -(Z-l) 

T]  "  U 


(2.1.1) 


F  (r),kx)  is  the  Coulomb  wavef  .action  (Abramowitz  and  Stegun  1972) 


-TIT] 


FX(T1'P)  =  A     S^  P"+Ie_ip  1F1(^l+in,2i+2,2ip), 

(2.1.2) 
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and  1F1  is  the  confluent  hypergeometric  function.   The  energy 
eigenvalues  of  H(l)  are  given  by 


H(l)|£>  »$!]£  \h    .  (2.X.3) 


With  these  functions  we  write  the  radiator  matrix  elements 
of  M(to)  (Eq.  1.6.7)  as 


i   \  3 

*<«>.     .    =-^-LfdteiAut 
nl.Q.mnln        2   J o 
i  1    L    L  h 


x    I  /   dfcj   /   d^2   <it1nJl1iai|vi(r,l)|^2nJl3m3> 

X3m3 

*<VV3I    .■*«1>"*TI<r.l).1«»t/|'.-»W»|^V.    (2.1.4) 


where  ng  is  the  perturbing  electron  density  and  \     is  the  thermal 
wavelength, 


*T 


■2Ttpti2a/2 


The  ket  |knlm>  is  the  product  of  perturber  and  radiator  kets, 
|k>|nlm>.   From  Eq.  (2.1.1)  we  see  that  the  k-angle  dependence  of  the 
perturber  wavefunction  <x|k>  can  be  written  explicitly  as 


<x|k>  =  I        I     Y*  (k)<x|kin>  .  (2.1.5) 

A=0  m=-X  X 


We  use  this  in  the  expression  for  M(cj)  ,  obtaining  four  £,m-sums  over 
perturber  angular  momenta.   The  k^-   and  k~-angle  integrals  each  are 
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integrals   of   two   spherical  harmonics  .     These   can   be   done   to   obtain 

Kronecker  deltas   which   reduce    the   four   x,m-sums    to   two,   which  we   label 

X.  ,m.    and   x. ,mr  .     We   then  have 
4     4  60 


M(o>)    ,  .  =  -  ^  T  dt   eiAa,t   f  dk,    J"  dk, 

nl.m,  ,ni„m„  _2  -1  o  ■'o        1   J  o        I 


nVl  •°*2,"2  f. 


where 

g  „      „    (k,,k_)  =  n  \Jc?k?    y      y      y 

n!    m      n       m    v     1  '     9'  s»    T    1     9  £  t,  L 


1    n  n        (k,  ,k.)   =  n  \  k,k„        y 

nJLra.  ,nJLra0v    1'    2J  e  fl   2      / 

i    I         2    2  X. 


~m„      JLm.      x,m, 
3   3        4   4        6   6 


x  <k1x4m4nx1m1|vi(r,l)|k2Jl6m6nx3ni3> 


x  <k2l6m6nx3m3|vi(r,l)|k1x4m4nl2m2>    .  (2.1.7) 


In  Appendix  B  we   reduce   G(k,,k2)    to   the  form 

4n  X e                                n-1  2n-2        » 

3  5  4      6 

(2x  +l)(2x.+l)(2x,+l)      JL  A     x,    -  x.    x.  A,    . 

„          3               4               6           r    1  3      5-i2  f    4      5      6 >  2 

(2x  +1) *>                   J  I  J 

V      5      ;                        0  0     0  0     0     0 


[£  dx  F     <„     k  *)F  „   (n2,k2x)A ,    „A  (x)]2    ,  (2.1.8) 

4  6  13   5 


where  X  . 

AxlV5(x>    =  C  dxr\2  ^x/V^  ~  4^)   V^V    ' 

X>  (2.1.9) 
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RnJ^Xr^  is  the  hydr°genic  radial  wavef unction  of  the  radiator,  and 
x<  (x>)  is  the  lesser  (greater)  of  xr  and  x,  the  radial  coordinates  of 
the  radiator  electron  and  perturbing  electron.   The  I      sum  is  over 
initial  radiator  angular  momenta,  the  I     sum  is  over  multipoles  of  the 
Coulomb  interaction,  and  the  Z^   and  X     sums  are  over  perturber  angular 
momenta.   The  symbol 


ml  m2  m3 
is  the  Wigner  3-j  symbol  (Edmonds  1957).   In  Appendix  C  we  express 

A|H  (x)  as  a  sum  which  is  convenient  for  numerical  evaluation. 
13  5 

The  time  integration  in  M(u>)  (Eq.  2.1.6)  can  be  done  by  adding  the 
convergence  factor  ie  (e>0)  to  Ago,  then  taking  the  limit  as  e-K): 

i«   ,.2       ,2. 

lim  f  dt   ei(A^fi£>t      e2m        *   '     2        - 
e*0     ° 


=   lim 


e+0     Aoo  +  ~  (k2  -  hh  +  ie 

=   '     Aaj+_^(k2.k2)     -  +  *«^  +i^l~  k2»>  (2.1.10) 

AtJ  +  2m   Ul        V 

where    P   stands    for    the    Cauchy   principal    part.      This    separates  M(u)    into 
real   and   imaginary   parts, 

1  co  co  e  2m      G(k,  ,k0) 

M   (Aa>)   =  —  P  /     dk,    /     dk 1 I t2   i    m 

R  «.2        Jo        1   Jo        2        .  -h        2  2        '  U.i.il) 


fi 


and  o*2,2 


AW  +  2m-(kl    "  V 

ph"k9 


1        G(k,,(k2  +  2jE^)l/2^ 


M   (Aa))   --«rdk     e        2m  i    i»    i   -       ft   * ^ 

1  fi3      °         *  ,k2  "  2mAco   1/2  '    ^>0' 

1  * 

(2.1.12) 
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M^Aco)  =  e  pt,'Awl  MjdAwl),  AuKO.  (2.1.13) 

The  function  MR(Aoo)  represents  a  dynamic  shift  of  the  spectral  line 
due  to  the  interaction  of  the  radiator  with  the  perturbing  electrons  of 
the  plasma,  M  (Aw)  represents  the  width  of  the  line  due  to  the 
perturbing  electrons.   We  use  the  fact  that  (O'Brien  and  Hooper  1974) 


2  2 

--^i    G(V(k2+^>1/2) 
P  J0  dk2  e  fi 2 2 =  °»    Au>0» 


Aw  +  ^  (kj  -  kj) 


to   write  MR(Aw)    in    the    form 


MR(Aco)    = 


2    2 

pfi   1c 


*2   Jo      lJo      2  -ft    .    9  ? 


* 9 o .    Ao)>0. 

A"  +  2m"  (k     "   k    ^ 

(2.1.14) 


This  integrand  is  not  divergent  so  the  principal  part  notation  is  not 
necessary.   Similarly,  for  Aw<0, 


mr(a.)  =  -L  #k  J> 

Tl 

Aw  +  —  (k2  -  k2) 
2m  k  1     2' 

(2.1.15) 


nume 
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Equations  (2  .1  .12)-(2  .1  .15)  for  M^Aw)  and  M  (Aw)  are  evaluated 

rically  and  the  results  are  used  in  Eqs .  (1.5.19)  and  (1.4.13)  to 

generate  line  shapes  .  We  will  present  the  results  of  these  calculations 

in  Chapter  III. 

II  .2   Angular  Momentum  Sums 

The  function  G(k1 ,k2) ,  which  appears  in  Eqs.  (2  .1  ,12)-(2  .1  .15)  for 

M(w),  contains  infinite  sums  over  the  perturber  angular  momenta,  JL  and 

4 

Jig .  Increasing  X   values  in  these  suras  correspond  to  increasing  separ- 
ation of  the  radiator  and  perturber,  so  for  I.    and  JL  greater  than  some 

4      6 

value,  say  V  ,  we  can  use  the  dipole  approximation  for  the  radiator - 
perturber  interaction  with  negligible  error  .  The  part  of  the  sum  in 
which  the  dipole  approximation  is  used  can  then  be  summed  exactly . 
We  separate  G(k^ ,k~)  into  two  parts, 


G(klSk2)  =  G(1)(k1,k2)  +  G(2\kvk2),  (2.2.1) 


where  G^   is  to  be  calculated  numerically  from  Eq  .  (2.1.8)  with  the 

upper  bound  V    on  the  l^   and  l^   sums,  and  G^  is  to  be  calculated  from 

a  dipole  approximation  of  Eq  .  (2.1.8)  for  the  remainder  of  the  I,    and 

4 

1&  suras.  The  dipole  approximation  for  Eq  .  (2.1.8)  is  obtained  by  re- 
stricting the  15  sum  to  the  dipole  term,  I     =  1 ,  and  replacing  Eq  . 
(2.1.9)  for  A^  z  x   (x)  by 

^r(n2_,y/2  ,2_  Ui+1)]l/26 


2Zx       l  Xl '*3+1         l  *i >X3  1 


(2.2.2) 
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We  obtain  this  last  expression  from  A      (x)  by  replacing  x 7x2  with 

2  13  5  <  > 

xr/x  ,  which  is  consistent  with  the  assumption  that  the  perturber  is  far 

from  the  radiator.  Since  j^-l  in  G(2),  we  can  use  the  equation  (Edmonds 
1957) 


l0  0  0j   l  i;    L(2JJ+3)(2A+1)J  (2.2.3) 


to  evaluate   the   3-j    symbols    in   G(2);    then  we   sum  over   I     to   obtain 

G(2)  (k     k   )   - 

n^ra^nJ^m^W    ~ 

,,34, 
4n   A  e        3na 

-^-har)    6i1i26m1m2(n2-^r1)fr(ki'k2)  <2'2^ 


where 


Wk2>   =  Jt+/{[/>FA-1(^.k1x)FA(r)2,k2x)  -L  ]2 


+   t/0   ^VV^WVV0  ij  ]2}    .  (2.2.5) 

x 

The   A-sum  in   this   equation   can  be   done   analytically   (Biedenharn   1956) 
with   the   result 

,,2    2 
2      2      2k1T11  -1 

fr(krk2)  =  r(k1+k2+-^i)[/odx  x"1  F1I-1(nlik1x)FJlt-1(T,2ik2«)]2 

_2rklk2    I1   +'HJ   j     I1   +"F   |    t/o  dx  x"1Fr(n1,k1x)Fr(n2,k2x)] 
x[/odxx"lFr-l(T|l'klx)Fi'-l{,,2'k2x)]     •  (2.2.6) 
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.(2) 


If  we  set  V    equal  to  one  in  G^  ^ ,  we  obtain  the  dipole  form  of 
M(w)  used  by  O'Brien  (1970)  and  Tighe  (1977).   O'Brien  showed  that 
f1(k1,k2)  is  related  to  the  free-free  Gaunt  factor,  g(k1,k2),  by 


VW^WW 


Karzas  and  Latter  (1961)  and  O'Brien  (1970,1971)  have  calculated  the 
free-free  Gaunt  factor  numerically .  They  reduce  the  integrals  in 
g(k^,k2)  to  sums  which  are  convenient  for  computer  evaluation.  The  work 
of  O'Brien,  which  was  for  singly  ionized  helium,  was  extended  by  Tighe 
to  radiators  of  arbitrary  charge  .  We  extend  their  work  to  evaluate 
fJl'(kl'k2)  f°r  V>1'     We  use  this  in  E(J«  (2.2.4)  to  obtain  G(2);  then  we 
calculate  M(w)  from  G  =  G^1^  +  G^  according  to  Eqs  .  (2  .1  .12)-(2  .1  .15)  . 
We  calculate  M(oj)  for  increasing  values  of  A'  to  determine  when  the 
dipole  approximation  in  g'  ^  gives  negligible  error,  i  .e  .  when  M(oo) 
becomes  independent  of  V  . 

II  .3  Electron  Correlations 


In  calculating  Eqs.  (2.1.12)  -  (2.1.15)  for  M].(Aoj)  and  M^(Aw),  we 
used  the  approximation  that  the  perturbing  electrons  do  not  interact 
with  each  other,  thereby  removing  effects  due  to  electron-electron 
correlations  .  Here  we  reintroduce  in  an  approximate  manner  the  effect 
of  these  correlations  on  the  line  shape  . 

These  correlations  have  been  shown  to  produce  a  screening  of  the 
radiator-perturbing  electron  interaction  (Hussey  et  al  .  1977).  This 
screening  is  most  significant  for  the  part  of  the  line  shape 
corresponding  to  times  which  are  long  compared  to  electron  relaxation 
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times,  that  is,  for  frequencies  inside  the  electron  plasma  frequency. 
Smith  (1968)  and  Hussey  et  al  .  (1977)  show  that  the  inclusion  of  corre- 
lations has  little  effect  on  M  (Aw)  for  |Aw|  >  ui   and  that  for 
IM  <  u   the  result  for  M^Aoj)   including  correlations  is  nearly  equal 
to  the  uncor related  result  evaluated  at  the  plasma  frequency,  M-.(w  ). 
Hence,  to  approximate  the  effects  of  correlations  on  fcL(Aw),  we  set 

M  (Aw)  =  MT(w  )  for  lAwl  <  w   . 
l        I  pe      '   '    pe 

To  approximately  include  electron  correlation  effects  in  >L,(Aw)  we 

replace  the  Coulomb  interaction  V^r.l)  in  G(k,,k2)  (Eq.  2.1.7)  by  the 

"X1/XD 
Debye-screened  interaction  V  (r,l)e      ;  then  use  the  screened 

G(k1,k2)  in  calculating  Eqs  .  (2.1.14)  and  (2.1.15)  for  M^(Au).  This 

approximate  treatment  of  correlations  follows  from  a  calculation  by 

Dufty  and  Boercker  (1976),  based  on  a  ring  approximation  treatment  of 

electron-electron  interactions  (Brout  and  Caruthers  1963),  in  which  they 

obtain  a  fully  screened  form  for  M(w) , 

M(w)  =  nTr1^(r,l;Aw)f(rl)[Aco  -  L(l)  -^(rl)]"1  x£(rl)f (r)"1 , 

(2.3.1) 

where  ^(r,l;Aw)  is  a  dynamically  screened  radiator-perturbing  electron 

interaction.  Then  they  proceed  to  show  that  for  Aw  <<  w   and  in  the 

pe 

non-degenerate  limit,  the  radiator-perturber  interactions  can  be  written 
in  the  Debye-screened  form 


2     -|x  -x  |/\     2  -x./Xn 
V'.l>--T-T-e   »  r   D-fe  »  D  .        (2.3.2) 

For  the  plasma  conditions  which  we  will  consider,  \     is  significantly 
larger  than  the  radiator  size  so  we  approximate  Eq .  (2.3.2)  by 
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Vg(r,l)  =  VjCr.l)  e  l      D    .  (2.3.3) 


Then,  using  the  approximations  of  Section  1.6  (except  for  the  neglect  of 
electron  correlations),  we  can  write  Eq .  (2.3.1)  as 


M(co)  =  nTr  V  (r,l)[Ao)  -  L(l)]  :V  (r,l)f(l)  .        (2.3.4) 

J.  o  S 


This  is  the  same  as  Eq .  (1.6.7)  except  for  the  screened  interactions. 
Following  the  steps  specified  in  Section  II  .1  we  can  obtain  from  Eq  . 
(2.3.4)  the  Equations  (2.1.14)  and  (2.1.15)  for  M-CAw),  except  that 
G(k1,k2)  will  contain  the  screened  interaction,  Eq .  (2.3.3). 

II  .4  Interaction  Matrix  Elements  Between  States  of 
Different  Principal  Quantum  Number 

So  far  we  have  assumed  that  matrix  elements  of  the  radiator- 

perturber  interactions  between  radiator  states  having  different 

principal  quantum  numbers  are  negligible  and  have  set  them  equal  to  zero 

(see  Appendix  A)  .  With  these  approximations  and  the  No  Lower  State 

Broadening  Approximation,  the  calculation  of  a  Lyman  line  shape  involves 

matrices  which  have  rows  and  columns  labeled  by  the  initial  radiator 

states  of  the  line  .  By  using  these  approximations  we  have  neglected 

quenching,  or  radiationless  transitions  between  radiator  states  of 

different  principal  quantum  number  caused  by  the  perturbing  electrons  . 

We  have  also  neglected  Stark  shifts  of  higher  then  first  order  in  the 

ion  microfield  and  the  changed  oscillator  strengths  due  to  mixing  of 

radiator  states  of  different  principal  quantum  number  by  the  ion 

microfield  . 
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For  a  given  line  we  expect  that  the  most  significant  of  the 
neglected  matrix  elements  are  those  between  the  initial  states  of  a 
given  spectral  line  (principal  quantum  number  =  n)  and  the  states  of  the 
nearest  adjacent  energy  level  (principal  quantum  number  =  n+1).   Here  we 
will  include  these  matrix  elements  in  the  calculation  of  the  n+1  Lyman 
line  shape.   The  matrices  involved  in  this  calculation  will  have  rows 
and  columns  labeled  by  radiator  states  of  principal  quantum  number  n  and 
n+1.   We  write  Eq.  (1.5.19)  for  J(u,e)  in  the  matrix  form 

J(o),e)=--Im  I   d   •  [to  -  (H(r)/-n-co  )  -  MCu)]?.1,  f,  ,  I. ,  t  ,    (2.4.1) 

11  11   1  AJ-  1  XI   1   xl 

where  H(r)  is  the  Hamiltonian  for  the  radiator  in  the  presence  of  the 
ion  microfield  e,    ho>  is  the  energy  of  the  radiator  ground  state,  and  i, 
i'  represent  radiator  states  having  principal  quantum  numbers  n  or  n+1. 
We  will  neglect  the  line  shift,  so  in  Eq.  (2.4.1)  we  have  set  B=0,  and 
we  will  set  M  (oj)  =  0. 

Following  the  derivation  of  Eqs.  (2.1.6)  and  (2.1.8)  from  Eq. 
(1.6.7)  we  can  write  M(w).  ,  as 

iAto   t 

*B)"lVl.Vrt"-7  n3V>e   ^   ^^ 

fl(k2_k2)t   _^i 
^^   3   2   6    ^  Gn1l1m1,n2V2,n3X3(kl'k2)      (2'4'2) 

where  Aw    =  oo-(w  -co  )  =  u-co   ,  fiw   is  the  unperturbed  energy  of  the 
states  with  principal  quantum  number  n,.   Here,  the  principal  quantum 
numbers  n^ ,  n2,  and  n3  can  have  the  values  n  and  n+1.   The  term  G(k1,k2) 
is  given  by 
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/      i    3   4 

4ne\     e  2n 


5  4      6 


(2V1)(2V1)(2X6+1)      *     X3  i         ^   Z5   I 

(2JS+1)       ""      lo    0    0  J  lo    0    0  j 


t  o  1    j    1   J   5 


r0  vvv)vvv)Vn,*1Wx)dx  •        (2-4,3) 

*  o  2   J    1    J   5 


The   term  A(x)    is   given  by 


,     5  6X     0 

V3W5W    =  ^  ^  S^(Xr)(-Vr  -  Hr"  K  *  <V    •      <2'4'*> 

X> 

We  use  Eq.  (2.1.10)  to  do  the  time  integration  in  Eq.  (2.4.2)  for  M(w), 

obtaining  for  the  imaginary  part  of  M(w): 


2   2 

2m 
e 


Ml(U)V,WA  =   "?     2,      Jodkl    Cdk2 
1112   2   2  ti     no^o 

x  ViV^vVa*1'*2'6^1"  h  (k2  -  ki>)  •      (2-4-5) 


We  rewrite  Aw    as  Aw  ,  =  Aw  -  w   ,  where  Am  is  the  frequency 
n_l      n_l         n.n        ^w  -1    j 

separation  from  the  unperturbed  line  center,  Aw  =  w-w   .   Then  we  do  one 

nl 

of  the  k-integrations  to  obtain 
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I         n^^^n^ 


2   2 

Pft  kl      g  (v      (\r2+  2mAw,  1/2- 

-«[Y    /-  dk    .'"ST"     ViVVzVW  :      x       * 


ft3    JL=i    °      x  fk2  +  ^  ]i/2 

3  l   1  -fi     J 


,    v      r00,,  2m      v   2       ft 

A3=0 


'1~1"V"2"2U12'"3'V 


V,»,.»,M,  fn,JL^k2+ Tf  (u,^»        K 


[}fy\w-t**)}m 


],  Aco>0;      (2.4.6) 


I         n^ij,^^^ 

pft2k2 
.n-1  -phlAwl    - 


-?n0W  2°    Sv,»2v2.v3(^+^-|,1/2'k2: 


3 


+  (second  term  of  Eq .  (2.4.6))],     AtXO;  (2.4.7) 

where  w1  =  u>     .  Here  we  have  broken  the  n^-sum  into  two  parts:  the 
first  term  in  the  brackets  of  Eqs  .  (2.4.6)  and  (2.4.7)  is  for  n3=n;  the 
second  is  for  n3=n+l  .  We  evaluate  Eqs.  (2.4.6)  and  (2.4.7)  numerically; 
then  use  these  results  in  Eqs.  (2.4.1)  and  (1.4.13)  to  obtain  line 
shapes.  Results  of  this  calculation  are  presented  in  Section  III  .2  . 


CHAPTER  III 
RESULTS 

III.l   Symmetric  Line  Profiles 
Here  we  compare  the  results  generated  using  this  full  Coulomb 
version  of  the  relaxation  theory  with  those  results  from  the  relaxation 
theory  in  which  the  dipole  approximation  is  made  for  the  radiator- 
perturbing  electron  interaction  (Tighe  1977;  Tighe  and  Hooper  1976, 
1978).   We  will  also  examine  the  relative  importance  of  various  multi- 
poles  of  the  Coulomb  interaction  in  M  (Ato)  by  restricting  the  A  -sum  in 
Eq.  (2.1.8)  to  the  multipoles  of  interest.   In  this  section  we  will 
neglect  the  line  shift,  which  is  given  by  the  terms  B  and  M  (Aw),  and 
consider  only  the  broadening  effects  of  the  electrons,  given  by  MfAw) 
(Eqs.  2.1.12  and  2.1.13).   We  will  examine  line  shifts  and  asymmetries 
in  section  III. 2. 

Comparison  of  Coulomb  and  Dipole  Results 

In  Figures  1  and  2  we  compare  Lyman  a  and  8  line  profiles  calcu- 
lated from  this  full  Coulomb  formalism  with  profiles  in  which  the  dipole 
approximation  was  used.   The  lines  appearing  in  these  figures  were 
calculated  for  a  plasma  of  Ar    ions  and  electrons  at  a  density  of  1024 
electrons  per  cubic  centimeter  and  a  temperature  of  800  eV,  conditions 
which  approximate  the  plasmas  of  some  recent  laser  implosion  experiments 
(Yaakobi  et  al.  1980).   For  the  Lyman  a   lines  presented  here,  the 
neglected  fine  structure  splitting  of  about  0.35  Ryd.  would  have  a 
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Figure  1   Comparison  of  an  Argon  Lyman  a  line  profile  calculated  from 
the  full  Coulomb  interaction  with  the  corresponding  profile 
calculated  from  the  dipole  approximation  for  that 
interaction  . 
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noticeable  effect  on  the  lines  (Lee  1979),  but  this  should  be  negligible 
in  the  case  of  the  p  lines.   Figure  1  shows  a  slightly  broader  central 
peak  for  the  a  line  calculated  from  the  Coulomb  interaction  than  for 
the  a  line  calculated  from  the  dipole  approximation.   In  Figure  2 
the  p  line  from  the  Coulomb  calculation  has  a  slightly  less  pronounced 
central  depression  than  does  the  corresponding  dipole  line.   The  Coulomb 
and  dipole  calculations  show  little  difference,  but  this  agreement  is 
partly  fortuitous.   Although  the  dipole  approximation  of  the  Coulomb 
interaction  neglects  broadening  due  to  the  other  multipoles  of  the 
interaction,  it  overestimates  the  dipole  (I  =D  part  of  the  Coulomb 
interaction  for  perturbers  which  are  less  than  a  few  times  the  average 
radiator  diameter  from  the  radiator.   This  nearly  compensates  for  the 
multipoles  neglected  in  that  approximation.   Figures  3  and  4  show  that 
Lyman  a   and  p  lines  calculated  from  the  dipole  approximation  are 
significantly  broader  than  those  calculated  from  the  dipole  part  of  the 
Coulomb  interaction. 

Although  we  have  removed  the  long  range  monopole  from  V,(r,l)  (Eq. 
1.6.2),  the  multipole  expansion  of  V^r.l)  still  contains  a  monopole 
contribution  from  the  part  of  the  perturber  wavefunction  which 
corresponds  to  penetration  of  the  radiator,  as  can  be  seen  from  Eq. 
(C2).   If  we  compare  a  and  p  lines  calculated  from  the  full  Coulomb 
interaction  with  lines  calculated  using  the  monopole  plus  dipole  part  of 
the  Coulomb  interaction,  the  lines  from  the  full  Coulomb  calculation  are 
only  slightly  broader,  indicating  that  for  these  conditions  the  most 
significant  broadening  of  these  lines  is  due  to  the  monopole  and  dipole 
parts  of  the  interaction. 


Figure  3   Comparison  of  an  Argon  Lyman  a  line  profile  calculated  from 
the  dipole  part  of  the  full  Coulomb  interaction  with  the 
corresponding  profile  calculated  from  the  dipole 
approximation  for  the  full  Coulomb  interaction . 
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Comparison  with  Impact  Theory 

In  Figures  5  and  6  we  compare  line  shapes  calculated  from  this  full 
Coulomb  relaxation  theory  to  line  shapes  calculated  from  a  full  Coulomb 
impact  theory  (Griem  et  al.  1979;  Kepple  1980).   This  impact  calculation 
approximately  accounts  for  the  screening  of  electron  fields  and  the 
finite  duration  of  perturbing  electron  collisions  through  an  appropriate 
choice  of  cutoffs  of  the  electron  impact  parameter.   Griem  et  al.  also 
include,  by  an  approximate  method  based  on  partial  wave  scattering  cross 
sections,  the  quantum  mechanical  effects  of  strong  electron-radiator 
collisions.   Figures  5  and  6  show  that  our  a  and  0  lines  are  slightly 
broader  than  those  of  Griem  et  al.   In  Figure  7  we  compare  our  result 
for  a  p  line  at  a  density  of  8.2  x  10  3  electrons/cm3  and  a  temperature 
of  800  eV  to  the  results  of  Griem  et  al.  at  1024  electrons/cm3  and  800 
eV.   We  obtain  close  agreement  on  the  line  wings,  the  part  of  the  line 
which  we  believe  to  be  the  most  reliable  for  experimental  density 
determination.   Hence,  in  this  range  of  plasma  conditions  our 
calculation  would  indicate  a  plasma  density  about  18%  lower  than  would 
the  calculation  of  Griem  et  al. 

III. 2   Shifts  and  Asymmetries 
So  far  we  have  neglected  sources  of  asymmetry  in  the  line 
profiles.   In  this  section  we  will  examine  two  significant  sources  of 
line  asymmetry:  the  line  shift  due  to  perturbing  electrons  and  radiator- 
perturber  interaction  matrix  elements  between  states  of  different 
principal  quantum  number. 


Figure  5   Comparison  of  an  Argon  Lyman  a  line  profile  calculated  from 

this  full  Coulomb  formalism  with  an  a  line  profile  calculated 
from  the  full  Coulomb  impact  formalism  of  Griem,  Blaha,  and 
Kepple  . 
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Line  Shifts 

Our  objective  here  is  to  examine  the  line  asymmetry  which  is  caused 
by  the  static  and  dynamic  shifts,  B  and  M  (Aw),  and  to  determine  the 
significance  of  the  dynamic  shift  in  comparison  to  the  static  shift.   We 
obtain  M^(Aw)  by  the  numerical  evaluation  of  Eqs.  (2.14)  and  (2.15). 
Several  calculations  of  static  shifts  due  to  mean  field  effects  of  the 
plasma  on  the  radiator  have  been  done  (Skupsky  1980;  Davis  and  Blaha 
1981;  and  references  cited  therein).   Although  some  uncertainties  remain 
in  these  calculations,  for  our  purposes  we  will  use  the  results  of 
Skupsky  (1980,  1982)  for  the  static  shift.   These  results  were  obtained 
by  the  numerical  solution  of  the  Schrodinger  equation  for  a  hydrogenic 
ion  in  a  self  consistent  potential  determined  by  the  nonlinear  Poisson 
equation  for  electrons  and  neighboring  ions. 

The  Lyman  B  lines  in  Figure  8  include  the  static  shift  (solid  line) 
and  the  static  plus  the  dynamic  shift  (dashed  line).   The  static  shift 
causes  an  asymmetry,  the  lower  energy  (red)  wing  having  the  greater 
intensity.   The  dynamic  shift  causes  a  slight  red  shift  of  the  line  but 
has  no  significant  effect  on  its  shape.   This  full  Coulomb  result  for 
the  dynamic  shift  is  significantly  smaller  than  that  obtained  in  the 
dipole  approximation  (Woltz  et  al.  1982). 

Interaction  Matrix  Elements 


In  Figure  9  we  compare  Lyman  6  lines  calculated  with  (solid  line) 
and  without  (dashed  line)  radiator-perturbing  Ion  and  radiator- 
perturbing  electron  interaction  matrix  elements  between  states  of 
principal  quantum  number  3  and  4  (see  Section  II. 4).   The  inclusion  of 
these  matrix  elements  causes  the  line  to  be  asymmetric,  with  the  blue 
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peak  more  intense  than  the  red  peak  and  the  red  wing  more  intense  than 
the  blue  wing.   To  distinguish  between  ion  and  electron  effects,  we  set 
the  ion  perturbation  matrix  elements  equal  to  zero  while  retaining  the 
electron  perturbation  matrix  elements,  obtaining  the  solid  line  in 
Figure  10.   This  line  is  symmetric,  and  it  is  slightly  broader  due  to 
quenching  than  the  dashed  line,  in  which  quenching  has  been  neglected. 
The  asymmetry  of  the  (3  line  in  Figure  9  is  primarily  due  to  quadratic 
Stark  shifting  of  radiator  energy  levels  by  the  ion  microfield.   To  shov; 
this,  we  have  included  the  quadratic  Stark  effect  in  an  otherwise 
symmetric  Lyman  p  line  calculation  by  adding  the  quadratic  shift  (Bethe 
and  Salpeter  1977)  , 


AE  =  "16  ao  (f-)V[l7nx-  3q'-  9m'  +  19] 

r 


to  the  diagonal  elements  of  the  matrix 

<3qm|Ao)  -  hL(Aw)  |  3q'm'> 

(see  Eq.  1.5.19).  Here  the  kets  |nqm>  represent  states  of  the  radiator 
calculated  in  parabolic  coordinates,  and  the  quantum  number  q  is  equal 
to  the  difference  of  the  quantum  numbers  n,  and  n~  used  by  Bethe  and 
Salpeter.   This  calculation  gives  a  Lyman  p  line  which  has  an  asymmetry 
nearly  identical  to  that  of  the  p  line  in  Figure  9,  with  the  only 
noticeable  difference  being  that  this  calculation  gives  a  slightly 
larger  peak  asymmetry.  The  calculation  of  Figure  9  does  not  give  the 
exact  quadratic  Stark  effect  since  it  does  not  include  radiator-ion 
perturbation  matrix  elements  between  radiator  states  of  all  principal 
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quantum  numbers ,  but  the  near  agreement  of  these  two  calculations  shows 
that  the  matrix  elements  between  states  of  n=3  and  4  are  responsible  for 
a  major  part  of  the  quadratic  Stark  shift.  The  asymmetry  of  Figure  9  is 
similar  to  the  asymmetry  obtained  by  Griem  (1954)  from  a  Hydrogen  Balmer 
(3  line  calculation  which  included  the  quadratic  Stark,  effect  . 


CHAPTER  IV 
CONCLUSION 

In  this  work  we  have  extended  the  line  shape  formalism  of  Tighe 
(1977)  and  Tighe  and  Hooper  (1976,  1978)  by  eliminating  the  dipole 
approximation  which  they  used  for  the  radiator-perturbing  electron 
interaction,  retaining  instead  the  full  Coulomb  radiator-perturbing 
electron  interaction.  We  calculate  the  line  width  and  shift  operator, 
M(u)),  to  second  order  in  this  interaction.  A  second  order  calculation 
of  M(o))  will  give  good  results  for  the  part  of  the  line  profiles  in 
which  we  are  interested,  i.e.  separations  from  line  center  less  than  or 
on  the  order  of  the  electron  plasma  frequency,  since  the  electron 
broadening  of  this  part  of  the  line  is  primarily  due  to  weak  electron- 
radiator  collisions  .  Our  full  Coulomb  calculation  gives  Lyman  a  and  8 
lines  which  are  only  slightly  broader  than  those  calculated  from  the 
dipole  approximation,  but  the  near  agreement  of  the  two  calculations  is 
partly  fortuitous  .  In  the  dipole  approximation,  broadening  due  to  the 
other  multipoles  of  the  interaction  is  neglected,  but  this  is  partially 
compensated  by  the  fact  that  the  dipole  approximation  overestimates  the 
dipole  part  of  the  Coulomb  interaction  for  perturbers  which  are  less 
than  a  few  times  the  average  radiator  diameter  from  the  radiator  . 

The  Lyman  a  and  8  line  profiles  which  we  have  obtained  from  this 
full  Coulomb  calculation  are  somewhat  broader  than  corresponding  line 
profiles  calculated  by  Griem  et  al .  (1979)  from  a  classical  path  impact 
theory  with  quantum  mechanical  corrections  for  strong  collisions  between 
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the  radiator  and  perturbing  electrons .  For  plasma  conditions  on  the 
order  of  10   electrons/cm  and  800  eV,  density  diagnostics  based  on  our 
Lyman  8  line  would  indicate  a  plasma  density  about  18%  lower  than  would 
those  based  on  the  8  line  of  Griem  et  al  . 

We  have  calculated  the  dynamic  line  shift  due  to  perturbing 
electrons,  obtaining  a  small  red  shift  of  about  0.03  Ryd  .  for  a  Lyman  8 
line  calculated  for  the  plasma  conditions  of  10   electrons/cm  and  800 
eV .  This  shift  is  significantly  smaller  than  that  obtained  in  a  similar 
calculation  in  which  the  dipole  approximation  for  the  radiator- 
perturbing  electron  interaction  was  used  .  The  effect  of  the  dynamic 
shift  calculated  here  is  considerably  smaller  than  that  of  the  static 
shift  (Skupsky  1980,  1982),  which  causes  an  asymmetry  in  the  8  line, 
with  the  red  wing  having  the  greater  intensity.  This  asymmetry  results 
from  different  static  shifts  of  the  various  angular  momentum  states 
which  contribute  to  the  line  . 

We  have  included  in  the  calculation  of  a  Lyman  8  line  the  ion- 
radiator  and  electron-radiator  perturbation  matrix  elements  between  the 
initial  states  of  the  radiator  ,* principal  quantum  number  =  3,  and  the 
nearest  adjacent  energy  level,  principal  quantum  number  =4.  The 
inclusion  of  these  matrix  elements  gives  an  asymmetric  line,  with  the 
blue  peak  more  intense  than  the  red  peak  and  the  red  wing  more  intense 
than  the  blue  wing  .  The  asymmetry  is  due  to  the  matrix  elements  of  the 
ion-radiator  perturbation,  which  cause  nonlinear  Stark  shifts  of  the 
radiator  energy  levels  .  The  electron-radiator  perturbation  causes  a 
slight  additional  broadening  of  the  line  due  to  quenching  . 

Sholin  (1969)  has  shown  that  when  the  quadratic  Stark  effect  causes 
a  significant  line  asymmetry,  the  quadrupole  interaction  between  the 
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radiator  and  the  ions  will  also  have  a  significant  effect  on  the  line 
profile.   Calculations  of  Balmer  and  Paschen  p  lines  which  include  the 
quadrupole  interaction  (Pittman  and  Kelleher  1981)  give  asymmetric  line 
profiles  with  the  blue  peak  being  more  intense  than  the  red  peak.   The 
quadrupole  interaction  should  cause  a  similar  asymmetry  in  the  Lyman  p 
line,  although  the  magnitude  of  the  asymmetry  for  the  spectral  lines  and 
plasma  conditions  of  interest  here  remains  to  be  determined.   This  and 
the  quadratic  Stark  effect  may  be  significant  causes  of  increased 
intensity  of  the  blue  peak  of  Argon  Lyman  p,  an  asymmetry  which  has  been 
inferred  from  some  recent  laser  implosion  experiments  (Yaakobi  1982). 


APPENDIX  A 
MATRIX  ELEMENTS  OF  M(w) 

Here  we  will  use  the  No-Quenching  Approximation  (NQA)  and  the  No 
Lower  State  Broadening  Approximation  (NLBA)  to  write  M(oo)  (Eq  .  1  .6  .4)  in 
the  matrix  form  given  by  Eq .  (1.6.7). 

We  see  from  Eq  .  (1.5.19)  that  since  M(co)  appears  in  the  inverse 
operator  [w-L(r)-B-M(<jo)  ]   ,  which  can  be  considered  as  a  power  series 
expansion,  it  acts  on 

[L(r)+B+M(cj)]kf(r)o\  k>0  . 

Here,  the  operator  f(r)  weights  initial  radiator  states  and  the 
operator  d  determines  transition  probabilities  between  initial  and 
final  radiator  states;  so  the  only  matrix  elements  of  f(r)d  which 
contribute  significantly  to  a  given  line  are  of  the  type  [f(r)d]  f, 
where  i  and  f  represent  initial  and  final  radiator  states  for  that  line 
(Griem  1974)  .  We  can  set  all  other  matrix  elements  of  f(r)d  equal  to 
zero  .  Let  us  define  the  operator 


QK    J    =  [L(r)+B+M(o))]  f(r)d.  (Al) 


All  matrix  elements  of  D    which  are  not  of  the  form  TV  '  are  set  equal 

if 

to  zero  .  We  will  show  that  as  a  consequence  of  this  matrix  structure 
of  D    ,  the  matrix  elements  of  M(cj)D    ,  when  restricted  by  the  NQA  and 
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the  NLBA,  can  be  written  as  (Smith  1966) 


[M(u)D(0)]if  =  I  MCw)^,  »l°,\    ,  (A2) 

->(k) 
where  M(co)  .  .  ,  is  given  by  Eq.  (1.6.7).   Next  we  will  show  that  D    has 

the  sane  matrix  structure  as  D    ,  so  Eq.  (A2)  holds  for  all  D    and 

the  matrix  form  of  M(w)  (Eq.  1.6.7)  can  be  used  to  calculate 

[w-L(r)-B-M(oj)]_1. 

If  we  let  u  and  v  represent  arbitrary  radiator  states,  the  matrix 

elements  of  M(to)D    ,  where  M(w)  is  given  by  Eq.  (1.6.4),  can  be  written 


iM(»)i$«\v--;W>°ia,t 

Tl 

-i(E  f+H(l))t/1i  ,n.   i(Ef+H(l))t/1i 

_2  T    -i(.t«o»«/,  f(1)B»,       /C.*^ 

4*-  ui  if   fv 

,  r„  ,„/,xi  ,,        ,rts   ifE,+H(l))t/"fi 
*•_  ui     if  fv 

-i(E,+H(l))t/fi       ...      i(E  .+H(l))t/fi 

+  I   5  .  e    i         f(l)D  "V  ,e   V  V.J,    (A3) 

.2  ,   ui  'if   fv'  v'vJ 

ifv'   r 

where  we  have  written  the  term  [o)-L(r)-L(  1)  ]    in  M(w)  as  a  time  trans- 
form and  we  have  replaced  the  Liouville  operators  by  the  corresponding 

commutators.   Here,  V   ,  represents  matrix  elements  of  V,(r,l). 
uu'   r  1 

Now  we  use  the  NQA  and  the  NLBA  to  simplify  this  expression.   In 
the  NQA,  we  assume  that  the  perturbing  electrons  do  not  cause  the 
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radiator  to  make  nonradiative  transitions  between  states  of  different 
principal  quantum  number;  that  is,  the  matrix  elements  V   ,  are  nonzero 
only  if  the  states  u  and  \xf   have  the  same  principal  quantum  number  . 
This  restricts  the  u's  to  ^  initial  states  and  the  v's  to  be  final 
states  in  Eq .  (A3)  .  We  will  examine  the  significance  of  quenching  in 
Section  2.4.  In  the  NLBA,  we  assume  that  there  is  no  broadening  of  the 
final  radiator  state  by  the  perturbing  electrons,  so  we  set  the  matrix 
elements  Vfff  equal  to  zero.  This  is  a  good  approximation  for  Lyman 
lines  since  the  final  (ground)  state  and  the  plasma  have  no  dipole 
interaction,  which  is  the  most  significant  multipole  in  broadening. 
With  these  approximations,  the  second,  third,  and  fourth  terms  in  brack- 
ets in  Eq .  (A3)  are  equal  to  zero,  and  the  first  term  can  be  written  as 


[M(u))3(0) 


if 


ifTr.J>   I  e1A"tVil..e-1«"t/\..1,f(l)5<?>e1H<1>t/,i       (A4> 
•ft  i'i" 


where 


Aw  =  w-(El„-Ef)/1i  . 


If  the  final  states  are  degenerate,  Eq .  (A4)  can  be  written  as  the 
matrix  product 

*(0)i     v  „/..n   *(0) 


[M(U)5^]lf  -  X  M(M)llf^f 


where 
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M(w)    .  _J^Tr   J-dt  £  eiAut     -iH(l)t/n      iH(l)t/* 

ii'     ,2   1  Jo    j\,      li  i  i'         v  ' 

fi  i 

(A5) 
We  must  have  degenerate  final  states  (or,  as  in  the  case  of  Lyman  lines, 
one  final  state)  to  go  from  Eq  .  (A4)  to  Eq  .  (A5)  since  the  Aio  in 
M(cj)   ,  contains  E^  . 

We  see  from  Eq  .  (A4)  that  the  only  nonzero  matrix  elements  of 
M(w)D    are  [M(w)D    ].f.   This  is  a  consequence  of  the  NQA,  the  NLBA, 
and  the  fact  that  the  only  nonzero  matrix  elements  of  D    are  D  .,  . 
Similarly,  it  can  be  shown  that  the  only  nonzero  matrix  elements  of 

[L(r)  +  B]D^°'  are  {[L(r)  +  B]ir0^}.f,  so  D(1)  only  has  nonzero  ele- 
ments D  f    .     By  induction,  the  only  nonzero  elements  of  D    for  all 
k  >  0  are  D   ;  so  Eq  .  (A4)  holds  for  D    as  well  as  D    .  Hence,  we 
can  use  Eq .  (A5)  to  calculate  M(w)  in  [w-L(r)-B-M(u)] 


APPENDIX  B 
EVALUATION  OF  G(k  ,k  ) 

In  this  appendix  we  reduce  the  function  G  .     .   (k,  ,k„) 

nJLra    ,nJ(.  m„      1      2 

(Eq.   2.1.7)    to  a  form  which  is    tractable   for  numerical   calculation.     We 
have 


3,  2,  2 


;nim     nim  <kl»k2)    =   Wlk2      I  I  I 


,    , .        1  m 
3  3       4  4     16 

x  (kimnla    |v   (r  ,1)  |kj  m  nlji  ><k  J  m  nLm    |V   (r,l)  |k  J,  i.nLa  >    . 

(Bl) 
The  expression  for  V^r.l),  Eq .  (1.6.3),  can  be  expanded  in  spherical 
harmonics  (Jackson  1975)  to  obtain 

2       2 
V  (r,l)  =-^ S_ 

lx"xrl    lxl 

X 
-  I  x  5    6A  0 

-  4-2    *     X     tiitiT  Hbr  -  -f^W   ■  (B2) 

Jc=0  m=-x  5 

X> 

where  x  and  x  are  the  positions  of  the  perturbing  electron  and  the 

radiator  electron,  and  x^  (xj  refers  to  the  lesser  (greater)  of  x  and 

x  .  With  this  expansion  and  the  expression  for  the  perturber 

wavef unctions ,  Eq  .  (2.1.1),  the  first  matrix  element  of  V1(r,l)  in  Eq . 

(Bl)  can  be  written  as 


<k1A.m.nA,m1  |V.(r,l)  |k„JLm.nJLm.>  = 
144   111       266   33 
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A  -A   ^VV  -aU4,k  )) 

8e  i      e  ) 


X  m     5 


x  (/ dx  nAwYV5w\m/^(Kiiffl/;r^v/r^m3^r)) 


(£  dx  x  (kjx)   (k2x)  ^  (Tij.k^F^  (n2,k2x)Ax     (x))     (B3) 
'4         6         13  5 


where  A,     e 

x  5     6A  0 


and  R  (x  )  is  the  hydrogenic  wavef unction  of  the  radiator.  In  Appendix 
C  we  obtain  a  form  of  Eq.  (B4)  convenient  for  computation.  The  integral 
of  three  spherical  harmonics  is  given  by  (Edmonds  1957) 


/  dx  Y    (x)Y.   (x)Y.   (x) 
Vl    h*2  *3m3 

<2A1+l)(2A.+l)(2A.l+l)  ,,„   JL  A-  A.   A,  A.  A, 

=  [— " £ — ]1/2  (  *   2   3H  l      2     3)  (B5) 

m  m  m«   0  0  0 


where 


■  h    *2  A3- 


is  the  Wigner  3-j  symbol.   We  use  this  equation  and  the  fact  that 


Y*  (x)  =  (-l)m  Y.  m(x) 
xm  x-m 


to  obtain 


2       A  -A       i(o(A    ,k  )-a(A    ,k, ) 
<k1A4m4nA1m1|V1(r,l)|k2A6m6nA3m3>   =^-i6      4e 
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m,+m,+m,. 
14      5 


*  I      (-D    l  5[(2JL+l)(2Jl„+l)(2JL+l)(2VI-l)]1/2 

x_ra_ 

x    /  Xl    *3   *5w   *4   *5   V/l    *3   VA   *5    A6-j 

-m.   m_  mc      -m.-m.   m,      0     0     0       0     0      0 
13     5  4     4     6 

x   Q  dx   ^   (n1,k1x)F      (,2,k2x)A  (x).  (B6) 

4  6  13   5 


Similarly,  we  obtain  for  the  second  matrix  element  of  V,(r,l)  appearing 
in  Eq.  (Bl) : 


2        JL-JL      i(a(x     k   )-a(x   ,k  )) 
<k2V6ni3m3|V1(r,l)|k1i4m4nI2m2>=^-i4     6e  "      ' 

m_+m,+m7  . 

x   X      (-1)  D      /[(213+l)(2Jl2+l)(2Jl6+l)(2Jl4+l)]i/Z 

x    ,h   l2   *7w   h   h   *4w*3   *2   *7v*6   *7    V 

-m_  m„  m.,      -n.-ii7  m,      0     0     0        0     0      0 
3      2      7  6      7      4 

x    /;  dx   F     (Vk2x)F      (VklX)A  (x).  (B7) 

6  4  3   2   7 


We  use   Eqs.    (B6)    and    (B7)    to   write   Eq.    (Bl)    as 


C  (k     k    1   = 

nx  m1,nA2m2v  1'  2* 

4tlp^Te  1/? 

2  I  I        [(2Jl1+l)(2A2+l)]i/Z(2Jl3+l)(2Jl4+l)(2Jl6+l) 

3   4   6        5   7 

X    Uo   dX   Fl(Vklx)FL(Vk2x)A«Ll(x)) 
4  6  13   5 


x    (/;  dx    F      (t,2.k  X)P      (T,     k  X)A  (X) 

6  4  3    2   7 
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/l   *3  *5w*4  A5  A6uX2  A3  ^^4  X6  A7> 
000        000        000        000 

m,  +m„+m ,  +ra ,  +m,.+m-. 
x      I  I      (-1)    13     4     6      5      7 

m3m4m6  n,5m7 

^  (  a2  ji3  i5^  *4  *6  x6^  a3  i2  A7^  ae  z7  x4^  ^ 

-m,  m_  mr   -m, -mr  m,.   -ra„  m„  m_   -m^-m^  m. 
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The  m-sums  in  this  equation  can  be  done  with  the  aid  of  the  following 
properties  of  the  3-j  symbols  (Edmonds  1957): 


(1)    JL  ,  JL  and  J!,  must  satisfy  the  triangle  inequality, 

\i  -SL    |  <  A  <  (A+JL),  and  mj+n^+m-j  must  equal  zero  in  order  for 
the  3-j  symbol 


>h  h  z3^ 


to  be  nonzero . 


(2)  The  completeness  relation: 

1  O  4.  1  O         *5  i 

£()(,)=   (2A.+1)        60    .,6        ,A(A.,A,,,A,)    , 
L      vm,    m.   m.  Mm.    m„   m' ;  3  A,AI   m.m*         1'    2'    3 

mm        123123  3333 

where  A(A  ,A  ,A  )  is  one  if  JL.JL.JL  satisfy  the  triangle 
inequality  and  zero  if  they  do  not. 

(3)  An  odd  permutation  of  the  columns  of  a  3-j  symbol  is  equivalent 

Jl  +A2+A3 
to  multiplication  by  (-1)        : 
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A  l2   *3>  m   (_*l+V"'S/2  ll   *3^ 

ml  m2  m3  m2  ml  m3 

(4)    Changing  the  sign  of  all  the  m's  is  equivalent  to  multiplication 
I  +1+1 
by  (-1)  l     l     J: 


■h    h   l3^  m   (_1)A1+A2+A3/-  h    h   l3- 


A  special  case  of  this  rule  is  that  if  m,=m2=mo=0,  A  +JL+JL-  must 
be  even  or  the  3-j  sumbol  will  be  equal  to  zero. 


First  we  do  the  m-  and  m,-  sums, 


m,+m,  X.    Xr   X.        X,    Xn   X,  mc     X.    X,    Ac  X,    X,    A, 

I     (-1)  4  6(   4  5  6)(   6  7  4}  =  (_1}  5  j(  4  6  5)(  4  6   7) 


=  (-1)  5(2It.+l)6    6^      A(*  A  A  )  .  (B9) 


Next  we  use  the  Kronecker  delta  in  m,-  and  -m-,  to  do  the  sum  over  m^ , 
then  the  m0  and  mc  sums  are 


m,+m,  X.    A.  I       A_  A,  A. 
I   (-1)  3   5(   '   3   5)(  2   3  5) 


=  (-D  X(-l)  L  3  5  I   (  *  3  5)(  2  3  5) 
m_+m,.  m.-m--ia_  m  -m_-m,- 

m  +X  +X   +A 
=  (-1)  3(2A.+1)  \  ,  6m   A(A.,AVAS)  .        (BIO) 

1       12   12 


Using  these  results,  we  write  Eq.  (B8)  as 
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4n  A.  e 
nl  e    ,ni  n       1*    2  2  JLJL      mm 

n-1         2n-2        »  (2JL+1)(2JL+1)(2JL+1)      JL    JL    JL      I.     JL    JL    _ 

-    I       I       Z  3   m*n — ^—t1   3   5)(4   5   6)2 

JL=0     JL«0     JL.JL-0  v**5T*'  0     0     0       0     0     0 

3  5  4     6 

x    [J~  dx   F^   (n^k^Fjj   (ri2,k2x)AJl  A  ^   (x)]2    .  (Bll) 

4  6  13   5 


This   is   the   desired   result,   Eq.    (2.1.8).      We   use    this   equation   (with  Eq. 
(C.8)    for   A(x))    to   calculate   G(k, .kj)    numerically,    then   the   results   for 
G(k,,k2)   are   used   in  Eqs.    (2. 1 .  12)-(2.1 .  15)    to   obtain   M(io). 


APPENDIX  C 
A  COMPUTATIONAL  FORM  FOR  A      (x) 

12  3 
Here  we  reduce  the  term  A      (x)  to  a  form  which  is  practical  for 

X1X2  3 
computer  evaluation.  We  start  with  Eq  .  (2.1.9), 

\V3(X)  =  C  dXr\  ^V^I^T  ~  HR  R»jl  <*r>«      (C1) 


where  x>.  (xO  is  the  lesser  (greater)  of  x  and  x  ,  which  are  the  radial 
coordinates  of  the  perturber  and  radiator,  respectively.  The  function 
R  (x  )  is  the  hydrogenic  radial  wavef unction  for  a  radiator  with 
nuclear  charge  Z  (Messiah  1961), 


_   ,   .    r   (n-A-l)l  ,1/2    r/2  I   _  2JH-1   ,   ,  ,__. 

Rni(xr}  =  * V  e      Xr  Ln-x-l  (xr}'  (C2) 

2n[(n+x)!]J  r  n  x  i   r 

na 
where  x  and  x  have  been  scaled,  in  terms  of     ,  and  L  is  the  Laguerre 

polynomial , 


l2*+1  (x)  =nT1(-nk       i(*+*-y-i3*k  (C3) 

u-MlX;   kfQ  k    i;   (n-x-l-k)!(2x+l+k)!k!  *  tCJ; 


If  we  substitute  Eqs  .  (C2)  and  (C3)  into  Eq .  (CI)  we  obtain 


I 

n-JL-1  n-x„-l  k,+k„+x,+x„+2     -x       x„  xo,0 

X  I  Cvu     fdxx      L      2      X      2       e     r(_4_--J_)        (C4) 

kf-0       k,-0  klk2      °      r   r  ^      V1  X      J 

12  x. 

69  > 


70 
where 


v2  ■ 


k,+k.  [  (n-A.-l ) !  (n-A.-l )  !  (n+JL  )  !  (n+A„) ! 


(-1)  1   2  _ 1     2   '    ~l'^"  -2" 

'  2n(n-l1-k1-l)!(n-I2-k2-l)!(2Jl1+l+k1)!(2A2+l+k2)!k1!k2!  *  K^> 


We  see  from  Eq.  (2.1.8)  that  A      (x)  will  be  multiplied  by  a  3-j 

1*2  3 
symbol  containing  I   ,    A  ,  and  A  ,  so  A  ,A  ,  and  A  must  satisfy  the 

triangle  inequality,  |  JL  — J2.  |  <  A  <  (A.+A.)  .   Also,  kj  and  k2  are 

greater  than  or  equal  to  zero.   Therefore  the  exponent  of  x  in  Eq.  (C4) 

is  always  positive,  and  we  can  use  the  definite  integrals 


J*  dx  x*  e  *r  =  m!(l  -  e  X  \    fA ,  m>0  ,  (C6) 

j=o  J- 

and 

/"  dx  x"  e   r  =  m!  e"X  £  fr   ,  m>0  ,  (C7) 

j=o  J- 


to  write  Eq.  (C4)  as 


n-A  -1  n-A  -1 

AA  JLA  (X)  =   ^      *     Ck  k 
1*2  3       kj-0   k2=0     12 

A  ,0;  _   (k  +k  +A  +A  +A.+2)   j 

x  I— nr(WWV2),(l"a    A  fr 

3  j=0  J 

x  J 

A  (k  +k  +A  +A  -A  +1)  j 

+  x  J(k  +k  +A  +A  -A  +1)!  e  X   J1   Z   X   Z  J    ~ 

^  J        j=0  J' 

V°  , .  t„xi   -x  (^l+k2+W2)  x^ 


(k^k^+A^)!  ex  lL      L      l      L         *y  ]  .         (C8) 
j=0  J  • 


This  equation  is  used  in  the  numerical  calculation  of  A      (x)  . 

Xl*2  3 


APPENDIX  D 
NUMERICAL  METHODS 

To  obtain  numerical  results  for  the  line  profiles,  we  calculate 

M  (Aw)  and  MfAw)  numerically,  then  use  these  results  in  a  line  shape 

program  by  Tighe  (1977)  .  This  program  calculates  J(w,e)  by  inverting 

the  matrix 


[Aw  -  (H(r)/fi  -  w  )  -  B  -  M(w)] 


(Eq .  1.5.19),  then  does  the  ion  microfield  integration  (Eq.  1.4.13)  to 
obtain  I(w)  . 

We  calculate  M  ( Aw)  from  Eqs  .  (2.1.12)  and  (2.1.13),  using  Eq . 
(2.1.8)  for  G(k,  ^2)  .  The  Coulomb  wavef unctions  ,  which  appear  in  Eq  . 
(2.1.8),  are  calculated  by  a  continued  fraction  expansion  and  recursion 
relations  (Barnett  et  al  .  1974)  .  For  small  values  of  the  argument  kx 
the  Coulomb  wavefunctions  are  calculated  from  a  series  expansion 
(Abramowitz  and  Stegun  1972)  .  We  then  use  a  trapezoidal  rule 
integration  to  do  the  x-integral  in  G(k, ^2)  .  The  3-j  symbols  in 

G(k, ykj)    are  calculated  by  a  Fortran  function  given  by  Vidal  et  al  . 

2  2 
(1970).  We  make  the  change  of  variables  fjti  k  /2m  =  y  to  put  the  k^- 

integral  of  Eq .  (2.1.12)  into  a  form  which  can  be  integrated  by  Gauss- 

Laguerre  quadrature  . 

We  calculate  M^  Aw)  from  Eqs.  (2.1.14)  and  (2.1.15).  The  kj- 

integral  is  done  by  Gauss-Legendre  quadrature  for  a  mesh  of  k2  values , 

then  a  Simpson's  rule  integration  is  used  to  do  the  k2-integral  . 
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